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Einstein equations for several matter sources in homogeneous, isotropic metric are shown to reduce 
to a second order nonlinear ordinary differential equation. An analysis of its solutions is made in 
an important case. 

I. INTRODUCTION 

<n: 

Exact solutions of the Einstein equations are difficult to obtain due to their nonlinear nature. In this paper we 
^vq ■ show that the system of equations for homogeneous, isotropic cosmological models with a variety of matter sources 
reduce to particular cases of the ordinary differential equation 
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where /3, c and n are constants. 

The problem of a causal viscous fluid with the bulk viscosity coefficient £ proportional to p 1 ! 2 corresponds to n = 1, 
c = 0, y oc H in the truncated Extended Irreversible Thermodynamics theory,[l[ and n = —1/r, c = 0, y oc H~ r in 
the full theory, where the relation between temperature and energy density is assumed to be of the form T ~ /o r -0 
Also, the behavior near the singularity, when the relaxation term is much more important than the viscous term in 
the transport equation, corresponds ton = l, c = 0, y oc H for generic power-law relation £ = ap m .{3J] For a perfect 
1— 1. fluid source, with an equation of state p = (7 — l)p, and a cosmological constant A, we recover Eq.([T]) with n = 1, 
c = — 7 A provided we derive twice the 00 Einstein equation 
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We also find Eq.([T]) with n = 1, c oc A, y oc H for a fluid with equation of state 3p = —p — C/a 2 , p = Clog a/a 2 , C 
a constant. For two scalar fields, one free and the other self-interacting with a potential V(<j>) — Voexp(— A<p), the 
change y — a A I 2 in the Einstein equation leads to Eq.([T]) with n = —6/A 2 and c = 0, as can be seen in the paper 
"Two-scalar field cosmologies" appearing in this volume. For a time decaying cosmological "constant", A ~ —H 3 A 
with y oc H , the case with n = 1, c = also arises. [4] 

Thus, it turns out to be of great interest to analyze Eq. ([1]) from the mathematical point of view. Its general solution 
will be studied elsewhere, and we concentrate here on the families of real solutions of the case n = 1, c = 0. 



II. ANALYSIS OF THE SOLUTIONS FOR n = 1, c = 

Unless j3 = 1/9, Eq.([T]), for n = 1 and c = has only two point Lie symmetries and it is not equivalent to a second 
order linear equation under a point transformation. 5] So, we consider the nonlocal transformation 

z = y 2 , V = J ydt (3) 
which turns Eq.([T|) into the equation of a damped linear oscillator 

d 2 z dz 

+ +2f3z = (4) 
ar\ ar\ 

and we obtain the general solution of in a parametrized form (t(r)),y(rf)). The real solutions of (0} can be classified 
as follows: 
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a. /3 < 1/8 (strong damping). 

2(77) = Cexp(A±r/) (5) 

2(77) =2Ccxp(-77/2) cosh ((577/2 + 0) (6) 

21(77) = 2C*exp(-77/2)sinh (6r)/2 + 4>) (7) 

b. /3 = 1/8 (critical damping). 

z(77) = Cexp(-77/2) (8) 

2(77) = £7(77 + 0)exp(-77/2) (9) 

c. (3 > 1/8 (weak damping) 

z(r]) = 2Cexp(-r)/2)sm(5r)/2 + (p) (10) 

where \± are the roots of the characteristic polynomial, 5 = |1— 8/3 1 1 / 2 , and C, are arbitrary integration constants. 
Through the transformation ([3]), both Lie point symmetries of (H} have a simple equivalent: £(77) is defined up to 
an arbitrary integration constant to , and this reflects the invariance of (QJ under t — > t + to ■ Also, the invariance 
2 — > z/A 2 (A 7^ 0), of ((4|) is equivalent to the symmetry transformation t — > \A\t, y — > y/\A\. Besides, the permutation 
between the two branches of z leads to the discrete symmetry transformation y — > —y, t — > —t. 

Whenever 2(77) has a zero, extremum or inflexion point at 771 , y(t) has a zero, extremum or inflexion point at 
t\ = 4(771). Besides, it can be seen that y is finite at any zero point, and so y(t) is odd in a neighborhood of all finite 
zero points (see below). 

There are two groups of solutions 2(77): 

i) Those that never vanish, i.e., ([5]), ([6]), and (|5J), so that we may choose 2(77) > for all 77. For these solutions y(t) is 
nonvanishing, and it is obtained from any of the two branches of y/z (depending on sign y) . 

ii) Those that have (at least) one zero point, i.e., (0, ©, and (fTD) . The requirement that 2(77) > cannot be satisfied 
on both sides of the zero point by the same solution (with a given value of C). Therefore, 2(77) gives rise to two 
solutions y(t), one for each sign of C. Since y(t) is odd (see bellow), these solutions are obtained by joining at the 
zero point both branches of <Jz. 

From we see that 77(7;) is even for odd y(t) and has extrema at the zero points of y(t). Then, for the group 
(i), rj(t) is monotonic, that is, there is a one to one mapping between the real axis 77 and some interval of the axis 
t. However , for the group (ii), the pair of branches at each side of the zero point correspond to different mappings 
between 77 and t. Besides, each singularity of y(t) (where 77(7;) diverges logarithmically) marks a boundary for the 
mapping 77 — > t. For solutions ©, ((7]), JSJ), and (fTU|) . t(rj) can be expressed in terms of a hypergeometric function. 
Only for /3 = 1/9, or —1, t(rj) can be inverted in closed form. 

Due to the symmetries of ([TJ, if y(t) is a solution, Ay(AAt) is also a solution, where At = t — to ■ In particular, 
if y(At) satisfies y = and y ^ at t = to , — y(— At) is also a solution with the same initial data. However, as ([1]) 
satisfies a Lipschitz condition, given these initial data the solution is unique. Thus, we conclude that y(At) is odd. 
Further, it is easy to see that y(t) is analytic at to , so that there is an interval where its Taylor series converges. As 
y(t) must be odd also about any further zero point within the interval of convergence of the series, it comes out that 
there are only two possibilities for an interval with a zero point: 

a) The interval contains only one zero point. 

b) There are infinitely many equispaced zero points; that is, y(t) is an oscillatory periodic function (the radius of 
convergence is infinite). 

Solutions which exhibit behavior (a) occur only for (3 < 1/8, while the solutions for /3 > 1/8 have behavior (b). So, 
for (3 < 1/8, y(t) has either one or no zero point in the interval where it is bounded. 
The solutions ([5]) lead to the two one-parameter families of solutions for (3 < 1/8: 

y±(t) = a/At, a± = -2/A±, A+ ^ 0; y+ = K, \+ = (11) 

and we wish to investigate small departures from them. Let us consider first the case when exp[(A_ — A+)?7 — 2(f>] <C 1. 
As A + > A_ , this occurs for any if 77 is big enough. Then to first order we get the approximated solution 



v(t) = l^( 1 +l Atr ), A+^0; y(t) = K(l + 7 eM-Kt)), A+ = 



(12) 
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where 7 oc exp(— 2(f)), r — 4 — a is the Kowalevski exponent and a — a + in this case.Q In the opposite case, that is 
for r\ — > —00, we get also (fl"2"b ). but now a — a_ and 7 cx exp(2</>). Whenever y(t) has a singularity, z(rf) diverges, 
and this occurs for rj — > 00 as well as for i] — > —00 if ft < 0. Hence, using (j3|), we find that any singularity is located 
at a finite time, and in a neighborhood of it, y(t) has the asymptotic form (fT2"b ) with r > 0. When /3 > 0, z(rj) — > 
and \t(ri)\ — > 00 for 77 — > 00. Then, y(t) vanishes at infinity with the asymptotic behavior (|12b). where a — a+ and 
r < 0. 

The two-parameter families of solutions arise from ©, ©, ©, and (JTUJ). We classify them in two groups: those 
which have a singularity at a finite time, and those which are regular for all time; and we give the main features of 
their behavior, 
a. Singular solutions. 
< (3 < 1/8, C > 0. 

al. They have a singularity, where the leading behavior is (jilt -), decrease monotonically and vanish at infinity with 
leading behavior (jllN -) unless f3 = 0, when they have a nonvanishing limit. 
13 < 0,C>0. 

a2. They have two singularities, one with leading behavior (fTD-) and the other with leading behavior (fTTl f ). There 
is a minimum (maximum) between them. 

a3. They have a zero point between two singularities, where the leading behavior is They increase monotonically 

and have three inflexion points. 
l3< 1/8, C<0. 

a4. They have a zero point between two singularities, where the leading behavior is (fTTM for (3 < 1/8 and 



y(t) - ± 



1 A-lnln(l/|At|) 



21n|At| (21n|Ai|) 
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for j3 — 1/8. They decrease monotonically. 
b. Regular solutions. 
< 13 < 1/8, C > 0. 

bl. They have a zero point between two extrema. They have three inflexion points and they vanish at infinity, with 
leading behavior (|11N-) for f3 < 1/8 and like a4., with the replacement ln(l/|At|) — > In | At| , for (3 — 1/8. 
f3 = 0,C>0. 

b2. They have a zero point at to and increase monotonically with a nonvanishing limit at infinity. 
•3 > 1/8. 

b3. They are oscillatory periodic, and its period and amplitude have a relation of the form AT 2 = f(5). The period 
diverges as (3 1/8+ and has the limit T 2 7 r 3 / 2 /(v / C|r(3/4)| 2 ), (3 — !• 00 . These results agree with the phase space 
analysis, and confirm the numerical simulations. Q 



III. CONCLUSIONS 



We have obtained the general solution of Eq.([T]) for n = 1 and c = in a parametrized form by means of Eq.Q 
The solutions have moving singularities and depending on whether these points are real or not, two groups of real 
solutions arise: the singular and the regular ones.[9j Both one-parameter families of solutions for (3 < 1/8 are singular, 
unless (3 — 0, when one of them turns into a constant. They coalesce for (3 = 1/8, and there is no real one-parameter 
family of solutions for f3 > 1/8. These one-parameter solutions give the leading behavior of the solutions about a 
singularity. Only for f3 = 1/9,0,-1 two-parameter solutions are functions on the complex plane and real solutions 
can be expressed in closed form. 

In general, the problem of the construction of explicit solutions of a given integrable nonlinear differential equation 
remains open. One direction along which one can attempt to proceed is linearization, i.e. the reduction of the equation 
to a linear ordinary differential equation, which is, by definition, integrable. Only for j3 = 1/9 Eq.([l]) possesses eight 
symmetries and is linearizable by a point transformation. On the other hand, the transformation @ linearizes it for 
any value of (3. Thus, although it has only two Lie point symmetries, it possesses eight nonlocal symmetries. We 
think that it is of utmost importance to study this kind of linearizing transformations, which have received up to now 
little attention. 
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